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Summary 

Direct  numerical  simulations  (DNS)  were  carried  out  of  bypass  transition  in  favorable, 
zero,  and  adverse  pressure  gradient  boundary  layers  in  the  presence  of  free-stream  turbu¬ 
lence.  A  time  sequence  from  the  ZPG  simulations  was  stored  in  order  to  obtain  a  number 
of  spots  for  the  study  of  the  calming  effect.  In  the  last  six  months,  additional  spots  have 
been  obtained,  in  order  to  increase  the  number  of  available  events  to  forty  two  (42)  spots. 
Various  flow  visualization  techniques  were  used  to  analyse  the  evolution  of  the  spots.  It 
was  observed  that  the  location  of  spot  inception  is  difficult  to  identify.  In  the  previous 
report,  we  analyzed  the  instability  of  the  streaky  boundary  layer.  We  herein  report  the 
use  of  the  Lyapunov  Exponent  in  order  to  determine  the  Lagrangian  Coherent  Structures 
(LCS),  which  form  the  spot.  We  also  report  on  our  effort  to  determine  the  average  spot 
inception  location,  spreading  rate,  and  speed  of  propagation. 


1.  Introduction 

1.1.  Bypass  transition 

Boundary  layer  breakdown  to  turbulence  is  generally  classified  into  two  distinct  types: 
The  first,  natural  or  orderly  transition,  is  prevalent  in  low  disturbance  environments  and 
proceeds  via  the  amplification  and  breakdown  of  Tollmien-Schlichting  (TS)  waves.  The 
second  type  is  the  ‘bypass’  scenario.  It  is  preceded  by  the  amplification  of  boundary 
layer  perturbation  streaks,  or  Klebanoff  modes.  The  streaks  can  amplify  to  levels  of  10 
-  15%  of  the  free-stream  velocity,  LGo,  and  breakdown  at  low  Reynolds  numbers,  hence 
‘bypassing’  the  natural  route.  The  bypass  scenario  prevails  in  conditions  of  free  stream 
turbulence  above  1%  of  the  free  stream  velocity,  and  is  therefore  important  for  many 
engineering  flows.  An  integral  part  of  the  bypass  process  is  the  secondary  instability  of 
boundary  layer  streaks  (Swearingen  &  Blackwelder  1987;  Andersson  et  al.  2001),  which 
leads  on  to  the  formation  of  turbulent  spots  (see  figure  1.1).  These  spots  are  sporadic 
patches  of  turbulence  which  grow  and  merge  to  sustain  the  upstream  edge  of  the  fully 
turbulent  boundary  layer. 

The  significance  of  turbulent  spots  is  widely  acknowledged,  and  they  display  a  variety 
of  interesting  phenomena.  For  instance,  the  area  trailing  the  turbulent  spots,  known  as 
the  calmed  region,  has  been  shown  to  have  a  stabilizing  effect  on  the  flow  with  significant 
implications  for  separation  delay  (see  for  example  Johnson  2001). 

1.2.  Influence  of  turbulent  spots  on  boundary  layer  separation 

In  adverse  pressure  gradient,  boundary  layers  are  prone  to  separation  and,  as  a  result,  loss 
of  lift  and  higher  viscous  losses.  One  example  is  the  flow  in  turbo- machinery,  particularly 
in  the  compressor  passages,  due  to  the  adverse  pressure  condition. 

Figure  2  is  an  illustration  of  the  instantaneous  laminar  flow  around  a  compressor  blade, 
from  direct  numerical  simulations  (Zaki  et  al  2009,  2010).  Contours  of  instantaneous  vor- 
ticity  capture  the  separated  boundary  layers  on  both  sides  of  the  blade.  On  the  pressure 
surface,  the  separation  rolls  remain  adjacent  to  the  surface.  The  dynamics  on  the  suction 
surface  are  more  complex.  The  interaction  of  the  separated  Kelvin-Helmholtz  (K-H)  rolls 
and  the  newly  formed  boundary  layer  is  shown,  at  four  phases  of  the  shedding  cycle,  in 
figure  3  (a  —  d).  It  is  seen  that  the  Kelvin-Helmholtz  vortices  convect  slightly  away  from 
the  wall.  However,  they  maintain  a  separation  region  in  their  shadow,  within  the  newly 
formed  boundary  layer.  This  interaction  is  similar  to  that  reported  in  Jones  et  al.  (2008). 

Separation,  however,  can  be  delayed  or  bypassed  if  the  laminar  separation  point  is 
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Figure  1.  A  boundary  layer  undergoing  bypass  transition  due  to  free-stream  turbulence.  Con¬ 
tours  of  the  instantaneous  streamwise  velocity  perturbation  are  contrasted  at  two  planes:  one 
inside  the  boundary  layer  and  one  in  the  free  stream.  The  former  is  characterized  by  elongated 
Klebanoff  distortions  and  sporadic  bursts  of  turbulent  patches.  The  domain  is  duplicated  in  the 
periodic  spanwise  direction. 


preceded  by  boundary  layer  transition  to  turbulence.  Indeed,  the  simulations  by  Zaki 
et  al  (2009),  and  in  a  later  paper  by  Zaki  et  al.  (2010),  demonstrate  that  both  migrating 
wake  sand  free-stream  turbulence  can  trigger  bypass  transition  upstream  of  separation. 
If  the  transition  process  is  complete  upstream  of  separation,  the  boundary  layer  remains 
attached.  A  time  sequence  which  captures  this  behavior  is  shown  in  figure  4.  When  the 
turbulent  spots  form  upstream  of  separation,  they  spread  as  they  propagate  downstream 
and  cause  local  re- attachment.  This  process  is  intermittent:  local  re-attachement  events 
are  separated  by  relaxation  of  the  boundary  layer  towards  the  laminar  separated  state. 

In  order  to  demonstrate  the  calming  effect,  the  tangential  velocity  is  plotted  at  a  point 
within  the  separation  region  versus  time  in  figure  5.  The  figure  captures  the  influence 
of  the  convected  turbulent  spot:  Initially,  a  more  pronounced  separation  is  observed  at 
the  spot  leading  edge,  followed  by  local  re-attachment  of  the  boundary  layer.  The  re¬ 
attachment  is  abrupt,  and  persists  beyond  the  period  of  spot  passing.  The  boundary 
layer  indeed  separates  again  after  the  spot  has  traveled  downstream,  but  the  relaxation 
to  that  state  is  relatively  slow,  and  is  associated  with  the  calming  effect  of  the  spot. 

1.3.  The  anatomy  of  turbulent  spots 

Previous  experimental  literature  have  lead  to  three  main  conclusions  concerning  the 
universality  of  the  shape  and  kinematics  of  the  turbulent  patches.  Firstly,  the  shape  of 
turbulent  spots  resembles  an  arrowhead  pointing  in  the  positive  streamwise  direction. 


7 


FIGURE  2.  Snapshot  of  the  vorticity  held  in  the  disturbance- free  simulation  TO. 

(b) 


Spanwise  vorticity 

0  '  '  '5V  '  'lOO 


Spanwise  vorticity 

'  '  '5V  '  'lOO 


(c) 


Figure  3.  Four  phases  of  the  shedding  cycle  from  the  disturbance- free  simulation  TO.  Total 
period  is  T  —  0.24,  which  corresponds  to  non-dimensional  frequency,  F  =  10 qojv /U^  —  189 


Secondly,  the  tip  of  the  spot  leading  edge  is  elevated  and  overhangs  laminar  unperturbed 
fluid.  Finally,  the  shape  of  the  spot  appears  to  be  self-similar  as  it  evolves  following 
approximately  a  conical  similarity,  which  was  proposed  by  Cantwell  et  al.  (1978)  and  later 
used  by  several  studies  including  Vasudevan  et  al.  (2001)  and  the  more  recent  Sabatino 
&  Smith  (2008).  This  conical  similarity  is  based  on  an  observation  that  the  spot  spreads 
in  a  linear  fashion  in  the  streamwise  direction,  with  an  approximately  constant  spanwise 
spread  angle  of  about  20  degrees.  The  trailing  edge  (TE)  and  the  leading  edge  (LE) 
attain  constant  streamwise  propagation  speeds.  The  LE  moves  at  speed  approximately 
equal  to  90%  of  the  free-stream  velocity,  while  the  trailing  edge  propagation  speed  is  on 
the  order  of  50%  of  the  free-stream  speed.  The  difference  in  propagation  speeds  causes 
an  elongation  of  the  spot  as  it  travels  downstream. 

The  mechanism  for  the  spreading  of  turbulent  spot  has  also  been  the  subject  of  previous 
investigations.  The  laminar-turbulent  interface  evolves  as  the  spot  convects  downstream 
and  increases  in  extent.  Two  possible  entrainment  processes  have  been  put  forward:  The 
first  process  is  ‘nibbling’,  which  involves  the  propagation  of  the  laminar-turbulent  inter¬ 
face  by  means  of  small  scale,  viscous  diffusion  of  vorticity.  This  mechanism  is  dominant  in 
the  wall-normal  spreading  of  the  spot  (Cantwell  et  al.  1978).  However,  in  the  horizontal 
plane,  the  spread  of  the  spot  is  via  the  second,  more  complex  mechanism.  This  involves  a 
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Figure  4.  A  time  sequence  showing  the  inception  of  a  turbulent  spot  on  the  suction  surface 
of  the  compressor  blade.  An  iso- surface  showing  separation  is  also  shown.  The  spot  convects 
downstream  and  spreads  downwards  towards  the  wall.  It  causes  local  re- attachment. 


destabilization  of  the  laminar  region  outside  the  spot  due  to  the  turbulent  eddies  within 
the  turbulent  patch  (Riley  &  Gad-el  Hak  1985). 

The  structure  of  turbulent  spots  has  been  described  as  a  series  of  large-scale  spanwise 
vortices  (see  Wygnanski  et  al.  1976;  Cantwell  et  al.  1978).  However,  recent  flow  visualiza¬ 
tions  have  yielded  a  more  complex  structure,  including  hairpin  vortices  and  streamwise 
oriented  streaky  structures  (Sabatino  &  Smith  2008).  These  hairpins  contribute  to  the 
lateral  growth  of  the  spot,  also  influencing  the  entrainment  process  at  the  overhung  region 
of  the  leading  edge  (Krishnan  &  Sandham  2007). 

The  random  nature  of  spot  formation  has  caused  difficulty  in  experimentally  study¬ 
ing  their  evolution.  In  order  to  overcome  this  problem  and  control  the  location  of  spot 
inception,  researchers  have  employed  various  methods  of  triggering  spots.  For  instance, 
Cantwell  et  al.  (1978)  conducted  experiments  in  a  water  channel,  and  initiated  spots 
by  tapping  the  wall  of  the  channel.  Gad-el  Hak  et  al.  (1981)  carried  out  experiments  in 
which  turbulent  spots  were  triggered  by  means  of  a  solenoid  valve  ejecting  fluid  through 
the  flat  plate  into  the  boundary  layer.  The  turbulent  spots  in  the  majority  of  these  stud¬ 
ies  can  be  considered  as  artificially  triggered.  Naturally-triggered  spots,  on  the  other 
hand,  are  sporadic  and  difficult  to  capture.  Jacobs  &  Durbin  (2001)  conducted  DNS  of  a 
zero  pressure  gradient  boundary  layer,  subjected  to  free  stream  turbulence  (FST)  only. 
The  spots  that  were  observed  were  caused  by  the  interaction  of  lifted  low-speed  streaks 
and  free-stream  eddies  at  the  edge  of  the  boundary  layer.  We  herein  discuss  our  recent 
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Figure  5.  The  tangential  velocity  near  the  suction  surface  wall  is  plotted  versus  time.  The  time 
series  is  taken  at  a  point  downstream  of  the  laminar  separation  location.  The  figure  shows  the 
influence  of  the  passing  turbulent  spot,  and  the  calming  effect. 


Figure  6.  Schematic  structure  of  the  vortex 
loops  produced  in  early  stages  of  spot  inception 
from  McAuliffe  &  Yaras  (2008). 


Figure  7.  Primary  Hairpin  vortices  along  with 
subsequent  vortex  loops  in  the  early  stages  of  a 
turbulent  spot  using  the  A2  criterion,  fromKr- 
ishnan  &  Sandham  (2006). 


Primary  hairpin  structure 


progress  on  applying  visualization  technique  based  on  the  Lyapunov  exponent  in  order 
to  track  the  evolution  of  the  spots. 


2.  Methods,  Assumptions,  and  Procedures 

2.1.  Direct  numerical  simulations  of  bypass  transition 

Despite  the  moderate  Reynolds  number  in  transitional  flows,  the  computational  cost  of 
DNS  is  still  large.  For  instance,  the  inflow  plane  to  the  DNS  domain  must  be  place  up¬ 
stream  of  the  leading  edge  of  the  plate  in  order  to  allow  sufficient  time  for  the  synthesized 
inlet  turbulence  to  adjust  and  become  a  solution  to  the  Navier- Stokes  equations.  In  ad- 
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Figure  8.  Schematic  of  the  DNS  domain.  The  domain  boundary  is  indicated  by  the  dashed  line. 
The  curvature  of  the  top  boundary  ensure  a  constant  streamwise  pressure  gradient  parameter. 

dition,  the  leading-edge  region  must  be  fully  resolved,  which  significantly  increases  the 
grid  requirements. 

An  alternative  approach  which  reduces  these  constraints  was  proposed  by  Jacobs  & 
Durbin  (2001),  and  is  also  used  in  the  current  investigation.  A  schematic  of  their  DNS 
domain  is  shown  in  figure  8.  The  computational  domain  is  denoted  by  the  dashed  line, 
and  commences  at  Re  =  80  which  is  downstream  of  the  edge  of  the  plate  and  upstream 
of  the  transition  region.  This  choice  of  inflow  location  reduces  the  grid  requirements  for 
the  simulations. 

The  inlet  perturbation  is  synthesized  as  a  superposition  of  continuous  Orr-Sommerfeld 
and  Squire  modes.  The  continuous  modes  are  vortical  eigenfunctions  of  the  linearized 
perturbation  equations.  Th  shape  of  these  modes  resembles  sines  and  cosines  in  the  free- 
stream.  Therefore,  they  can  easily  be  used  to  synthesize  a  particular  turbulence  spectrum, 
with  randomized  phase.  Inside  the  boundary  layer,  the  modes  decay  as  they  approach 
the  wall.  Their  decay  is  consistent  with  the  solution  of  the  linear  perturbation  equation, 
and  hence  do  not  need  a  significant  downstream  length  within  the  computational  domain 
to  adjust  to  the  solution  of  the  Navier-Stokes  equations. 

The  domain  size,  in  terms  of  inflow  J99,  is  525  x  30  x  32  in  the  streamwise,  wall- 
normal,  and  spanwise  directions.  The  number  of  grid  points  in  each  of  the  coordinates 
was  1793  x  193  x  193,  respectively.  The  top  boundary  of  the  computational  domain  is 
curved.  The  curvature  is  designed  to  ensure  a  constant  pressure  gradient  parameter, 
A e  =  (62 /^(dUoo/dx),  or  equivalently  a  constant  Falkner-Skan  parameter,  /?,  where 
UOQ(x)  =  Four  such  DNS  were  carried  out,  at  A <9  =  {2,  0,  —  2,  —  3}  x  10-2. 

All  four  simulations  impose  a  zero  pressure  gradient  profile  at  the  inlet,  which  evolves 
downstream  into  the  correct  pressure  gradient  parameter. 

An  instance  of  the  boundary  layer  disturbance  field  is  shown  in  figure  9.  The  two  panes 
show  the  streamwise  u-  and  wall- normal  u-perturbations  inside  the  boundary  layer.  The 
R-perturbations  show  the  large-amplitude  streaks,  or  Klebanoff  distortions.  Both  panes 
show  a  turbulent  spot,  and  the  downstream  fully-turbulent  boundary  layer. 

2.2.  The  free- stream  turbulence 

The  literature  on  bypass  transition  has  clearly  demonstrated  the  importance  of  clearly 
documenting  the  free-stream  turbulence  conditions  when  comparing  the  influence  of  flow 
parameters,  for  instance  pressure  gradient,  on  transition  location.  In  our  DNS,  the  free- 
stream  turbulence  intensity  exhibits  the  same  time  history  among  all  simulations.  This 
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Figure  9.  A  boundary  layer  undergoing  bypass  transition  due  to  free-stream  turbulence.  Con¬ 
tours  of  the  instantaneous  streamwise  (top)  and  wall-normal  (bottom)  velocity  perturbations. 
The  u— perturbations  clearly  shows  the  elongated,  high-amplitude  Klebanoff  distortions.  Both 
panes  show  the  turbulent  spot,  and  downstream  fullty-turbulent  boundary  layer.  The  domain 
is  duplicated  in  the  periodic  spanwise  direction. 


ensures  that  spot  inception  and  spreading  rates  can  be  compared  in  a  meaningful  manner. 
The  turbulence  intensity  from  all  four  simulations  are  shown  in  figure  10.  It  is  clear  that, 
independent  of  the  pressure  gradient,  the  turbulence  in  the  free-stream  decays  similarly 
in  downstream.  Any  changes  in  the  transition  process  can  therefore  be  attributed  to  the 
mean  pressure  gradient,  rather  than  the  intensity  of  the  free-stream  forcing. 

The  average  location  of  transition  onset,  and  the  transition  length,  are  generally  ob¬ 
tained  from  the  time-averaged  skin- friction,  C/,  which  is  shown  in  figure  11.  All  four 
curves  start  at  the  zero  pressure  gradient  level,  since  the  inflow  is  a  Blasius  profile  plus 
free-stream  turbulence.  Downstream,  the  profiles  adjust  to  the  respective  skin  friction 
levels,  for  the  laminar  Falkner-Skan  solution.  Further  downstream,  transition  to  turbu¬ 
lence  takes  place.  The  figure  shows  that  adverse  pressure  gradient,  as  expected,  causes 
early  transition  to  turbulence,  and  a  shorter  transition  extent  in  comparison  to  the  zero 
and  finally  adverse  pressure  gradient. 

The  length  of  the  computational  domain  was  deemed  insufficient  for  the  zero  and 
favorable  pressure  gradient  simulations.  For  this,  and  other  reasons  discussed  below,  a 
new  simulations  was  carried  out  during  the  first  six  months  of  the  project.  The  new 
simulation  extends  the  computational  domain  to  x/5  =  600  instead  of  x/5  =  525,  and 
therefore  ensures  that  the  transition  process  is  complete  far  upstream  of  the  exit  plane. 

2.3.  Extended  simulation  domain 

The  time-averaged  skin  friction  curve  in  figure  11  suggest  the  streamwise  extent  of  the 
domain  is  insufficiently  long  for  the  zero  and  favorable  pressure  gradient  simulations.  In 
order  to  ensure  domain  size  independence,  the  streamwise  extent  of  the  domain  must 
be  extended  and  the  simulations  repeated  for  those  conditions.  It  was  therefore  decided 
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Figure  10.  Time- averaged  turbulent  intensity,  Tu,  for  various  values  of  A,  the  pressure 

gradient  parameter. 


Figure  11.  Time-averaged  skin  friction,  Cf ,  for  various  values  of  A,  the  pressure  gradient 

parameter. 

that  a  longer  domain  size  is  need  in  order  to  ensure  that  the  spreading  of  the  turbulent 
spots  is  unaffected  by  the  outflow  boundary  conditions. 

In  addition,  after  discussion  with  Dr.  John  Clark  at  the  AFOSR-JPR  in  Long  Beach, 
California  (2007),  it  was  agreed  that  the  number  of  turbulent  spots  from  the  previous 
simulation  might  be  limiting  in  the  convergence  of  the  proposed  statistical  study.  There¬ 
fore,  a  longer  time  period  of  the  DNS  simulation  must  be  stored  in  order  to  accumulate 
sufficient  turbulent  spots  for  converged  statistics. 

New  simulations  were  carried  out  for  the  zero  and  favorable  pressure  gradient  simu¬ 
lations.  The  computational  domain  for  these  simulations  extends  to  x/Sqq  =  600.  The 
time- averaged  skin  friction  curves  for  the  new  domain  are  shown  in  figure  12.  The  stream- 
wise  extent  is  now  sufficiently  long  to  ensure  independence  of  the  boundary  layer  evolu¬ 
tion  from  exit  boundary  conditions,  for  the  zero  pressure  gradient  simulation.  The  FPG 


13 


Figure  12.  Time-averaged  skin  friction,  C/,  in  the  zero  pressure  gradient.  This  new 
simulation  has  a  longer  streamwise  domain. 


simulation,  on  the  other  hand,  is  still  transitional  near  the  exit  plane.  As  a  result,  the 
late  stages  of  transition  are  not  reliable  since  they  can  be  affected  by  the  exit  boundary 
conditions. 

Snapshots  of  the  full  velocity  field  were  stored  throughout  the  simulation,  every  two 
convective  time  units.  The  total  number  of  snapshots  is  4000,  and  they  correspond  to 
the  flow  traveling  the  transition  length  more  than  36  times.  It  is  therefore  expected 
that  sufficient  spots  will  be  available  for  post-processing.  In  order  to  reduce  the  storage 
requirements,  the  DNS  fields  are  only  stored  every  two  time  units,  based  on  the  convective 
time-scale.  In  addition  only  every  third  point  of  the  DNS  field  is  included  in  the  output 
field.  The  requires  storage  from  the  zero  pressure  gradient  simulation  is  still  0.9  Tera- 
Bytes. 


2.4.  Visualization  of  the  turbulent  spots 

The  challenge  in  identifying  the  turbulent  spots  is  two- fold.  Firstly,  the  spots  are  spo¬ 
radic  in  nature;  their  inception  location  and  time  are  unknown  a  priori.  In  DNS,  the 
simulation  can  be  restarted  when  elevated  instantaneous  skin-friction  levels  are  detected 
at  the  wall,  and  the  spot  evolution  is  recorded.  The  second  difficulty  is  due  to  the  pres¬ 
ence  of  high-amplitude  disturbances  throughout  the  boundary  layer.  For  instance,  the 
amplification  of  boundary  layer  streaks  can  interfere  with  the  ability  to  detect  the  spot. 
This  is  demonstrated  in  figure  13  where  the  time  evolution  of  a  turbulent  spot  is  shown 
using  contours  of  ^—perturbation  velocity.  Instead,  use  of  vertical  perturbation  velocity 
contours,  shown  in  figure  14,  is  likely  to  result  in  better  spot  detection. 

Various  spot  detection  algorithms  have  been  tested  during  the  first  phase  of  the  project. 
Each  method  was  assessed  in  terms  of  its  ability  to  accurately  describe  the  shape  of  the 
spot,  its  spread  angle,  and  its  propagation  speed.  A  number  of  the  employed  visualization 
techniques  were  based  on  iso-surfaces  of  physical  flow  quantities.  These  methods  shared 
one  drawback,  which  relates  to  the  dependence  of  spot  size  on  the  threshold  of  the  selected 
iso-surface.  An  alternative  approach  based  on  the  Lyapunov  exponent,  and  which  does 
not  share  that  drawback,  is  also  considered  in  this  section. 
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Figure  13.  Instantaneous  contours  of  streamwise  velocity  fluctuation  (—0.2  <  u  <  0.2)  in  ZPG 
bypass  transition.  The  plane  is  located  at  y/So  =  0.74,  and  the  frames  translate  in  the  flow 
direction  at  half  the  free-stream  velocity. 


Figure  14.  Instantaneous  contours  of  wall-normal  velocity  fluctuation  (—0.05  <  v  <  0.05)  in 
ZPG.  The  plane  is  located  at  y/So  =  0.74,  and  the  frames  correspond  to  the  same  time  instants 
as  in  figure  13. 
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2.5.  Computation  of  the  Lyapunov  exponent 

Time-periodic  dynamical  systems  can  be  characterized  by  their  fixed  points  and  seper- 
atrices  which  distinguish  regions  with  different  dynamical  behavior.  For  example,  the 
seperatrix  on  the  phase  plane  of  the  simple  pendulum  system  seperates  the  stable  modes 
(pendulum’s  angular  position  always  less  than  7 r  rads)  and  the  unstable  modes  (pendu¬ 
lum  rotating  all  around  the  pivot  point).  However,  when  a  system  is  time  dependent, 
classical  methods  of  finding  these  invariant  manifolds,  formally  called  Lagrangian  Coher¬ 
ent  Structures  (LCS),  do  not  easily  apply  since  the  flow  map  of  the  system  changes  over 
time  (Haller  &  Yuan  2000). 

One  practical  way  of  identifying  LCS  in  time-dependent  dynamical  systems  is  by  us¬ 
ing  a  Finite  Time  Lyapunov  Exponent  (FTLE)  field.  In  qualitative  terms,  the  FTLE 
expresses  the  time  average  of  the  maximum  expansion  rate  of  a  pair  of  particles  advected 
in  the  flow  (Shadden  et  al  2005).  LCS  boundaries  can  be  identified  as  “ridges”  of  high- 
intensity  lines  in  the  FTLE  field.  The  method  was  explored  by  Haller  (see  Haller  &  Poje 
1998;  Haller  &  Yuan  2000)  and  revisited  by  Shadden  et  al  (2005).  The  main  advantage 
of  LCS  analysis  is  that  orderly  structures  can  be  revealed  for  a  finite  time  interval  for 
highly  chaotic  flows.  These  structures  would  not  easily  appear  through  traditional  tools 
of  analysis  (Shadden  et  al.  2005). 

LCS  analysis  has  been  applied  to  a  few  applications  in  fluid  dynamics.  A  popular  study 
was  related  to  contamination  accidents  which  can  potentially  affect  the  currents  close  to 
Monterey  Bay,  California  (Lekien  &  Leonard  2004).  Other  applications  include  studies 
of  2D  turbulent  structures  (Haller  &  Yuan  2000),  hairpin  vortices  (Green  et  al.  2007), 
vortices  of  bio-inspired  flows  (Shadden  et  al.  2005),  and  a  fully  turbulent  channel  flow 
(Green  et  al.  2007).  In  addition  Shadden  et  al.  (2005)  present  how  the  the  FTLE  field 
can  directly  reveal  the  separation  profile  of  the  flow  over  an  airfoil,  while  Lipinski  et  al. 
(2008)  have  presented  how  the  FTLE  field  can  be  used  to  study  2D  vortex-shedding  over 
an  airfoil. 

The  computational  implementation  for  obtaining  the  LCS  is  based  on  an  extensive 
treatise  by  Shadden  et  al.  (2005).  First,  let  D  C  Rn  be  the  domain  of  a  fluid  flow  with 
position  vector  x(t;to,Xo )  E  D  and  t  E  /  C  R.  A  dynamical  system  can  be  defined 
completely  by  the  following,  equations 


x(t]t0,x  0)  =  v(x(t]to,x0),t)  (2.1) 

^{to •>  to ^0^)  — ^0  (2.2) 

The  flow  map  can  be  defined  as 

4>t00+t(x)  =x(t;t0,x0),  (2.3) 

which  describes  the  trajectory  of  a  particle  starting  at  Xq  and  reaching  x  at  time  to  +t. 
Equations  2.1  and  2.2  are  solved  using  backward  integration  over  a  time  interval  T[nt. 
Backward  integration  is  used  to  identify  attracting  (unstable  for  time  independent  flows) 
manifolds,  while  forward  integration  can  be  used  for  repelling  (i.e  stable)  manifolds  over 
time  (Shadden  et  al.  2005).  For  the  flat-plate  boundary  layer,  backward  integration  is 
used  since  our  analysis  focuses  on  the  evolution  of  unstable  manifold  types. 

In  order  to  describe  how  neighbouring  particles  repel  or  attract  along  their  trajectories, 
we  introduce  the  matrix  M  (xij))  for  the  neighbourhood  of  a  particle  with  initial 
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Figure  15.  Relative  stretching  of  neighbouring  particle  trajectories.  Adapted  from  the 
discussion  in  Shadden  et  al.  (2005). 


position  a?o, 


dx 


(2.4) 


The  finite  time  version  of  the  Cauchy  deformation  tensor  for  the  same  particle  has  the 
form, 


A  |  i,j=MTM. 


(2.5) 


Maximum  stretching  of  the  neighboring  particles  corresponds  to  the  maximum  eigenvalue 
of  the  deformation  tensor.  This  eigenvalue  is  related  to  the  FTLE  of  particle  Xij  over 
Tint  according  to, 

°to+  (Xid)  =  77F  I  ln(\/^max)-  (2-6) 

Mint  | 

By  evaluating  a  for  all  points  Xij  in  the  domain,  an  intensity  FTLE  field  plot  can  be 
produced,  and  the  LCS  is  identified  as  lines  of  high  a. 


2.6.  Example  of  spot  visualization  using  the  Lyapunov  exponent 

The  time  evolution  of  a  turbulent  spot,  captured  using  the  Lyapunov  exponent,  is  shown 
in  figures  16  and  17.  The  sequence  spans  a  time  interval  AT  =  140  time  units,  based 
on  the  free- stream  mean  speed  and  the  inlet  boundary  layer  thickness.  The  top  view 
(figure  16)  shows  a  plane  at  distance,  y/So  =  1  from  the  wall,  where  So  is  the  boundary 
layer  thickness  at  the  inlet  of  the  computational  domain.  The  shown  region  of  the  flow 
domain  captures  the  inception  of  the  turbulent  spot  and  its  downstream  propagation. 
The  contours  of  the  Lyapunov  exponent  clearly  mark  the  turbulent  patch.  However,  some 
of  the  streaky  structures  within  the  boundary  layer  are  also  visible  in  the  figure.  Edge 
detection  schemes  are  still  required  in  order  to  identify  the  edges  of  the  turbulent  spots. 

In  figure  16,  the  leading  and  trailing  edges  of  the  spot  convect  at  different  speeds.  The 
former  is  moving  downstream  at  a  higher  speed,  and  hence  the  spot  grows  in  streamwise 
extent.  The  propagation  speed  of  the  spot  ‘center’  is  approximately  0.5  —  0.6£/oo,  but 
can  only  be  approximately  inferred  from  the  contour  plot.  In  addition  to  the  longitudinal 
spreading,  the  spot  also  spreads  in  the  spanwise  direction.  A  spread  angle  with  reference 
to  the  spot  trailing  edge  can  be  measured,  and  will  be  reported  in  subsequent  sections. 

Figure  17  shows  the  side  views  which  correspond  to  the  same  time  instants  as  in  figure 
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Figure  16.  Top  view  of  the  boundary  layer  showing  spot  evolution  using  contours  of  the 

Finite  time  Lyapunov  exponent. 
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X 

Figure  17.  Side  view  of  the  boundary  layer  showing  spot  evolution  using  contours  of  the 
Finite  time  Lyapunov  exponent.  The  snapshots  are  at  the  same  time  instants  as  in  figure  16. 
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16.  The  edge  of  the  boundary  layer  is  marked  in  the  figure,  and  is  defined  as  the  bound¬ 
ary  layer  95%  thickness.  The  side  views  clearly  demonstrate  that  the  spot  is  initiated 
near  the  free  stream,  due  to  an  outer  instability  of  the  streaky  base  flow.  This  is  the 
mode  of  breakdown  reported  by  Jacobs  &  Durbin  (2001).  Later,  Vaughan  &  Zaki  (2010) 
investigated  the  origin  of  the  outer  instability.  They  carried  out  linear  stability  analy¬ 
ses  of  streaky  boundary  layers,  when  the  base  streaks  are  unsteady.  They  demonstrated 
that  a  spanwise  sub-harmonic,  parametric  resonance  can  lead  to  fast  growing  waves.  The 
growth  rate  of  the  instability  was  shown  to  increase  with  the  unsteadiness  of  the  base 
flow,  which  is  the  case  in  the  current  simulations  because  the  Klebanoff  distortions  are 
generated  due  to  free-stream  vortical  forcing.  The  time  sequence  in  figure  17  also  shows 
that,  in  the  spot  region,  the  local  boundary  layer  thickness  is  increased  The  thickening 
of  the  boundary  layer  is  due  to  the  enhanced  momentum  mixing.  As  the  spot  convects 
downstream,  it  spreads  vertically  and  impinges  onto  the  wall.  This  impingement  is  fol¬ 
lowed  by  a  strong  enhancement  in  the  turbulent  intensity  of  the  spot,  and  faster  lateral 
spreading.  As  a  result,  the  spread  angle  of  the  spot  is  expected  to  depend  on  the  distance 
from  the  wall,  especially  in  the  early  period  of  spot  development. 

In  the  computation  of  the  Lyapunov  exponent,  two  parameters  must  be  prescribed, 
namely  the  integration  period  and  time-step.  Figure  18  shows  an  independence  study 
of  the  integration  time  step.  The  total  backward  integration  period  for  the  computation 
for  this  spot  was  fixed  at  80  convective  time  units,  based  on  the  free-stream  velocity 
and  the  inlet  J99.  Three  integration  time-steps  were  used,  At  =  {0.5,  1.0,  2.0}  time 
units.  The  Lyapunov  exponent  recomputed  at  the  same  y— location,  and  the  results  are 
in  agreement. 

Based  on  the  results  of  the  Lyapunov  exponent,  the  propagation  speed  of  a  turbulent 
spot  was  computed  since  inception.  Three  values  are  reported  in  figure  19  versus  time: 
the  speed  of  the  leading-edge,  mid-point,  and  trailing-edge  of  the  turbulent  spot.  The 
leading-edge  speed  is  approximately  90%  of  the  free-stream  velocity.  The  mid-point  of 
the  spot  convects  downstream  at  70%  of  the  free-stream  speed.  Finally,  the  trailing  edge 
of  the  spot  convects  at  approximately  50%  of  the  free-stream  velocity.  These  values, 
computed  at  y/S  =  0.04,  are  consistent  with  other  measures  in  the  literature  based  on 
the  wall  signature. 

In  addition,  the  spread  angle  of  the  turbulent  spot  can  be  computed  as  a  function 
of  time  since  inception.  The  results  are  shown  in  figure  20  for  two  planes,  one  near  the 
wall,  y/S  =  0.04,  and  the  second  higher  in  the  boundary  layer,  y/S  =  0.4.  The  spot 
width  is  larger  in  the  latter  case.  This  is  due  to  the  fact  that  the  spots  originate  near  the 
free  stream,  and  touch  down  onto  the  wall  later  in  time.  Therefore,  the  wall  signature 
is  initially  smaller  than  the  projection  of  the  full  spot  onto  the  horizontal  plane.  Later 
in  time,  the  spread  angle  in  both  the  near-wall  and  the  higher  plane  converge.  These 
results  stress  that  measurements  based  on  the  wall  signal,  for  instance  wall- stresses,  do 
not  convey  the  extent  of  the  spot,  particularly  at  the  early  stages  of  spot  spreading. 

The  above  results  were  based  on  the  analysis  of  one  particular  spot.  In  order  to  obtain 
reliable  measures  of  properties  such  as  propagation  speeds  and  spread  angles,  multiple 
spots  must  be  processed.  In  the  next  section,  a  time  map  of  the  spot  inception  is  pre¬ 
sented.  Using  edge  detection  algorithms,  spot  characteristics  will  be  obtained  from  the 
time  map.  In  addition,  the  conical  transformation  proposed  by  Cantwell  et  ai  (1978)  will 
be  introduced  and  used  in  the  following  sections  in  order  to  study  the  similarity  between 
the  various  patches  of  turbulence  recorded  in  our  direct  numerical  simulations. 
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Figure  18.  Effect  of  the  integration  time-step  on  the  Lyapunov  exponent. 
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Figure  19.  The  propagation  speed  of  a  turbulent  spot  is  plotted  versus  the  time  since  inception, 
o,  the  speed  of  the  leading  edge;  A,  the  speed  of  the  middle-point;  □  the  speed  of  the  trailing 
edge. 


3.  Results  and  Discussion 

3.1.  Analysis  of  turbulent  spots  from  DNS  database 

In  order  to  obtain  a  more  reliable  measure  for  the  average  speeds  and  spread  angles,  a 
similar  analysis  was  performed  for  all  the  spots  that  appear  in  a  time  period  of  7,940  time 
units.  The  recorded  42  spots  emerged  at  different  locations  and  times.  A  time  history 
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Figure  20.  The  spreading  angle  of  a  turbulent  spot  is  plotted  versus  the  time  since  inception, 
o,  the  spread  angle  at  y/S  =  0.4;  A,  the  spread  angle  at  y/S  =  0.04. 
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Figure  21.  Time  of  spot  inception  versus  spot  identifier. 


of  inception  is  shown  in  figure  21.  The  approximately  linear  rate  of  appearance  of  spots 
corresponds  to  one  spot  per  190  time  units,  based  on  Uqo  and  £99.  According  to  the 
time-averaged  skin  friction  profile  (see  figure  12),  the  length  of  the  transition  region  is 
approximately  100  length  units.  Based  on  a  propagation  speed  of  the  center  of  the  spot 
of  the  order  of  0.5E/ooj  the  lifetime  of  a  turbulent  patch  is  approximately  200  time  units. 
This  lifetime  is  consistent  with  the  rate  of  inception  (one  spot  every  190  units)  shown  in 
figure  21. 

It  is  important  to  note  that  the  spots  do  not  develop  in  exactly  the  same  pattern  as  the 
first  spot  shown  in  figures  16  and  17.  Instead,  they  differ  in  terms  of  inception  location 
and  interference  of  two  (or  more)  spots  developing  close  to  each  other.  A  histogram 
and  the  cumulative  distribution  function  of  spot  inception  with  respect  to  the  different 
streamwise  locations  are  given  in  figures  3.1a  and  6,  respectively.  Both  the  mean  and  the 
median  are  close  to  x  =  260,  measured  relative  to  the  inlet  of  the  computational  domain, 

%  inlet  —  33. 

The  time  map  of  the  spots  is  shown  in  figure  23.  It  is  a  contour  plot  of  the  Lyapunov 
exponent  at  a  particular  (y,z)  location,  plotted  as  a  function  of  (x,t).  The  results  in 
figure  23  are  recorded  at  (y,z)  =  (0.1,30).  This  line  along  the  x-coordinate  is  the  axis 
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Figure  22.  (a)  Histogram  showing  the  frequency  of  spot  inception  for  a  range  of  streamwise 
locations;  ( b )  The  empirical  cumulative  distribution  function  showing  the  percentage  of  the  total 
number  of  spots  that  were  observed  upstream  of  a  given  streamwise  location 


of  symmetry  of  some  of  the  turbulent  spots.  Since  spot  inception  is  a  function  of  space, 
similar  (x,  t)  maps  were  required  at  the  axes  of  symmetry  of  the  various  spots. 

A  zoomed-in  view  of  figure  23  is  show  in  figure  24(a).  The  latter  focuses  on  the  de¬ 
velopment  of  a  single  spot,  and  illustrates  the  linearity  of  growth  for  the  trailing  and 
leading  edges  of  the  patch  of  turbulence.  The  leading  and  trailing  edges,  denoted  LE  and 
TE  respectively,  are  identified  using  edge  detection  techniques.  The  horizontal  distance, 
xle  ~  %t e,  is  the  streamwise  length  of  the  spot  along  its  axis  of  symmetry.  In  addition, 
the  slopes  of  the  LE  and  TE  lines  are  the  propagation  speeds,  and  are  computed  using 
linear  regression.  The  propagation  speeds  for  the  marked  spot  are  shown  figure  24(b), 
where  the  coordinates  have  been  shifted  by  to  and  xo,  the  inception  time  and  location  of 
the  turbulent  spot. 

The  constant  speeds  of  propagation  of  the  leading  and  trailing  edges  (see  figure  24 b) 
suggest  that  the  spot  length  is  linearly  proportional  to  time.  Whether  all  spots  have  a 
constant  speed  of  propagation  is  examined  using  the  non-dimensional  variables, 


x 


* 


X  —  Xo 
Xf  -  Xo 


(3.1) 


t-tp 
tf  -  to 


(3.2) 


where  xo  and  to  corresponding  to  the  virtual  origin  while  x /  and  t /  correspond  to  the  last 
reliable  trace  for  a  given  spot.  Figure  25a  —  b  show  the  non-dimensional  coordinates  for 
few  of  the  turbulent  spots.  The  results  demonstrate  that  the  leading  and  trailing  edges 
indeed  propagate  at  constant  speeds  independent  of  the  choice  of  the  spot.  However,  the 
actual  speed  of  propagation  does  not  need  to  be  the  same  for  the  data  to  collapse  on  the 
line  of  unit  slope. 

The  above  methods  for  determining  the  virtual  origin  and  the  corresponding  prop¬ 
agation  speeds  were  applied  to  twelve  of  the  42  spots.  The  selected  turbulent  patches 
were  most  amenable  to  analysis  using  edge  detection  over  a  sufficient  time  period  of  their 
evolution.  The  remaining  spots  were  not  included  in  the  analysis  due  to  their  merging 
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Figure  23.  x-t  time  map  of  the  FTLE  domain  defined  by  y/5  —  0.1,  z  =  30  and  x/5  G  [80,  520] 


in  the  streamwise  or  spanwise  directions  with  adjacent  patches  of  turbulence,  or  due  to 
their  late  inception  which  resulted  in  early  merging  with  the  fully  turbulent  boundary 
layer  downstream.  The  relevant  results  are  given  in  table  1  and  are  in  good  agreement 
with  the  existing  literature,  apart  from  the  TE  speed  which  appears  to  be  slightly  lower 
than  half  Uoq  on  average. 


3.2.  The  conical  similarity 

Despite  the  random  nature  of  the  spots,  previous  work  has  sought  a  universal  description 
of  the  evolution  of  these  turbulent  patches.  In  particular,  Cantwell  et  al.  (1978)  proposed 
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Figure  24.  (a)  Time  map  of  a  particular  turbulent  spot,  showing  the  detection  of  the  boundaries 
which  define  the  leading  and  trailing  edges.  The  distance  xle  —  xte  corresponds  to  the  total 
streamwise  length  of  the  spot  along  its  axis  of  symmetry,  (b)  The  propagation  speed  of  the 
leading  and  trailing  edge  as  a  function  of  time. 
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Figure  25.  Testing  the  conical  approximation  using  x*  and  t*  for  (a)  the  leading  edge  and  (b) 

the  trailing  edge 


the  notion  of  conical  similarity.  The  conical  transformation  rescales  spots  at  different 
times  during  their  development  onto  the  coordinate  axes  defined  according  to, 


x  —  Xo 
Uoo(t  -  to) 


(3.3) 


Uoo(t  -  to) 
z  -  Zp 
Uoo(t  -  t0) 


(3.4) 


(3.5) 


26 


Spot  ID 

Inception 

Time 

UT  E /Uoo 

(±0.005) 

Uc/Uoc 

(±0.005) 

Uc/Uoo 

(±0.005) 

Ot-avrg 

(±0.05) 

3 

554 

0.86 

0.68 

0.49 

19.0° 

6 

1200 

0.91 

0.70 

0.49 

18.2° 

9 

1630 

0.86 

0.67 

0.48 

21.0° 

10 

1954 

0.83 

0.65 

0.47 

23.3° 

13 

2462 

0.88 

0.68 

0.47 

19.8° 

19 

3330 

0.87 

0.68 

0.48 

20.2° 

20 

3520 

0.89 

0.70 

0.50 

21.1° 

26 

4704 

0.82 

0.64 

0.46 

23.0° 

33 

5952 

0.93 

0.73 

0.52 

19.5° 

36 

6684 

0.88 

0.68 

0.47 

20.3° 

38 

7190 

0.88 

0.67 

0.46 

20.1° 

40 

7520 

0.89 

0.69 

0.48 

18.2° 

mean 

std.  dev. 

0.876 

0.030 

0.681 

0.020 

0.480 

0.020 

20.31° 

1.62° 

Table  1.  Statistics  of  spot  kinetics  based  on  a  smaple  of  12  spots  for  y/S  —  0.1 


The  new  coordinates  are  the  normalised  propagation  speeds  at  a  specific  location  (x,  t ) , 
measured  relative  to  the  reference  location  (xo,£o)-  The  reference  point  (xo,£o)  corre¬ 
spond  to  a  virtual  origin  of  an  assumed  linear  growth.  This  may  not  coincide  with  the 
actual  inception  time  depicted  in  the  x  —  t  map  since  linear  growth  does  not  appear 
instantaneously  at  inception.  In  the  paper  of  Cantwell  et  al  (1978),  the  transformation 
effectively  collapsed  the  patches  of  turbulence,  supporting  the  view  of  universality  of  spot 
kinematics,  and  self-similar  evolution.  Those  spots,  however,  were  generated  by  boundary 
forcing;  the  bottom  wall  of  their  channel  was  perturbed  in  order  to  induce  their  forma¬ 
tion.  In  the  current  simulations,  the  effectiveness  of  the  transformation  in  collapsing  the 
naturally-occurring  spots  should  be  addressed. 

By  definition,  £  for  the  leading  and  trailing  edges  should  be  constant  over  time,  ac¬ 
cording  to  the  conical  similarity.  The  value  of  £  should  be  equal  to  the  dimensionless 
propagation  speed  of  the  respective  edge.  Figure  26  demonstrates  that  the  lead  and  trail¬ 
ing  edge  speeds  of  the  various  spots  are  indeed  constant,  and  collapse  using  the  conical 
transformation. 

The  main  advantage  of  using  the  transformation  is  to  collapse  the  edges  of  all  the 
turbulent  patches,  at  all  instants  of  their  lifetime,  into  a  single  self-similar  shape.  In 
figures  27  and  28,  the  edges  of  various  turbulent  spots  are  plotted  in  the  (£,  £)  and  (£,  rj) 
planes,  respectively.  The  former  is  the  plane  view  and  the  latter  is  the  side  view  of  the 
spots,  in  conical  coordinates. 

In  the  plane  view,  approximate  agreement  in  the  shape  is  observed  (figure  27).  The 
mean  value  of  the  spread  angle  also  agrees  with  the  statistical  results  based  on  edge 
detection  in  the  (x,  t)  map  presented  earlier.  However,  the  agreement  of  the  spot  shapes 
in  the  conical  coordinates  is  not  sufficient  to  assert  that  a  universal  spot  shape  can  be 
identified.  The  same  applies  to  the  side  view,  in  the  (£,77)  plane  (figure  28).  The  figure 
captures  the  overhang  shape  better  than  the  region  of  the  spot  near  the  free  stream. 
Nevertheless,  the  transformation  gives  a  reliable  approximation  for  the  basic  kinematics 
of  turbulent  spots. 


27 


Figure  26.  The  conical  coordinate  £  as  the  dimensionless  propagation  speed.  Spot  10  05  Spot 
13  square;  Spot  19  o;  Spot  36  V5  Spot  40  A 


Figure  27.  x-z  edge  profile  of  three  spots  transformed  into  conical  coordinates.  Oimin  —  18.69° 
and  Oimax  —  21.93°  based  on  the  statistics  of  the  spread  angle  for  the  al  12  Spots 


3.3.  Spot  structure 

The  above  discussion  focused  on  the  outline  of  the  turbulent  patches.  Here,  the  vortical, 
turbulent  structures  within  the  spot  are  inspected.  Often,  the  turbulent  spot  is  considered 
to  consist  of  an  agglomeration  of  hairpin  vortices.  Vortex  identification  criteria  were 
applied  in  order  to  visualize  the  vortex  structures  within  the  patches  of  turbulence.  In 
particular,  the  A2  criterion  (Jeong  &  Hussain  1995)  was  used,  where  A2  is  the  second 
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Figure  28.  x-y  edge  profile  of  three  spots  transformed  into  conical  coordinates 
eigenvalue  of  S 2  +  fi2,  and 


S  =  i[Vu  +  (Vu)T] 

(3.6) 

n  =  f[Vu-(Vu)T], 

(3.7) 

The  vortical  structures  within  the  turbulent  spots  were  visualized  using  iso-surfaces  of 
the  Lyapunov  exponent  and  A2  =  —0.015.  A  comparison  of  their  top  views  is  shown  in 
figure  29.  The  general  spread  angle  of  the  turbulent  patch  is  captured  by  both  criteria. 
The  edges  of  the  spot  are  marked  by  a  closed  line,  which  highlights  the  overlap  between 
the  two  visualization  techniques.  This  region  of  the  spot  corresponds  to  the  main  vortices 
which  engulf,  or  entrain,  the  surrounding  fluid  into  the  turbulent  spot.  Away  from  the 
marked  regions,  and  in  particular  in  the  middle  of  the  spot,  there  is  a  lesser  degree  of 
agreement  between  the  two  criteria.  This  mismatch  is  more  evident  in  figure  30,  where  a 
three-dimensional  view  of  the  spot  is  shown. 

In  figure  30,  the  agreement  of  the  two  criteria  in  capturing  the  edges  of  the  spot  is  clear. 
In  the  middle,  the  iso-surface  of  A2  captures  the  presence  of  a  row  of  symmetric  hairpin 
vortices.  However,  the  finite  time  Lyapunov  exponent  does  not  captures  these  vortical 
structures.  The  disagreement  is  most  likely  due  to  the  choice  of  the  backward  integration 
time  in  the  computation  of  the  LCS,  which  was  sufficiently  long  in  order  to  capture  the 
overall  structure  of  the  spot.  As  such,  it  exceeds  the  lifetime  of  the  hairpin  vortices  which 
form  and  breakdown  on  a  shorter  time  scale.  Therefore,  at  least  two  backward  integration 
times  must  be  used  in  the  computation  of  the  finite-time  Lyapunov  exponent,  in  order  to 
capture  the  dynamics  which  contribute  to  the  evolution  of  the  spots:  The  first  integration 
time  should  correspond  to  the  global  evolution  of  the  spot,  and  the  second  should  address 
the  formation  and  breakdown  of  the  eddies,  or  vortices,  within  the  patch  of  turbulence. 

4.  Conclusions 

The  interaction  of  free-stream  turbulence  with  a  laminar  boundary  layer  leads  to  by¬ 
pass  transition  to  turbulence.  This  phenomenon  was  simulated  using  direct  computa- 
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Figure  29.  Spot  33:  top  view  showing  the  iso-surface  of  the  finite  time  Lyapunov  exponent  in 
green,  and  the  iso-surface  of  A2  =  —0.015  in  blue.  The  two  closed  lines  on  the  sides  mark  the 
cores  of  the  vortical  structures  where  the  Lyapunov  exponent  LCS  and  the  A2  isosurface  overlap. 


Figure  30.  Spot  33:  three-dimensional  view  showing  the  iso-surface  of  the  finite  time 
Lyapunov  exponent  in  green,  and  iso-surface  of  A2  =  —0.015  in  blue. 


tions  of  the  Navier-Stokes  equations.  The  numerical  simulations  capture  the  proceedings 
of  transition:  Klebanoff  streaks  amplify  within  the  boundary-layer  shear,  followed  by 
their  secondary  instability  and  finally  the  formation  of  turbulent  spots  which  spread  and 
merge  downstream.  A  time  sequence  spanning  7,  940  time  units  was  stored  during  the 
simulations.  This  period  captured  the  evolution  of  42  turbulent  spots.  Twelve  of  these 
turbulent  patches  were  clearly  identifiable,  without  interaction  with  neighboring  spots. 
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In  order  to  study  the  evolution  of  the  turbulent  spots,  various  visualization  techniques 
were  considered,  for  example  based  on  the  three  components  of  perturbation  velocity 
and  vorticity.  The  streamwise  velocity,  and  also  the  wall-normal  vorticity  are  not  well 
suited  for  inspection  of  the  spot  shape  due  to  the  high-amplitude  streaks.  Instead,  the 
other  components  of  the  perturbation  field  can  be  used.  One  major  drawback  of  this 
approach,  however,  is  the  subjective  definition  of  the  spot.  In  the  literature,  a  level  of 
perturbation  velocity  is  often  used  to  define  the  spot,  for  instance  ||Y||  >  3%.  It  was 
demonstrated  that  the  spot  spread  angle  is  sensitive  to  this  definition  by  comparing 
three  values  of  the  threshold.  This  measure  must  also  be  adjusted  in  the  vicinity  of  the 
wall  since  the  perturbation  velocities  decays  in  the  near-wall  region.  The  use  of  the  finite 
time  Lyapunov  exponent  in  order  to  determine  the  structure  of  the  turbulent  spots  was 
investigated.  It  was  demonstrated  that  the  method  can  capture  the  turbulent  patches, 
without  any  clutter  due  to  the  Klebanoff  streaks. 

An  x  —  t  map  of  the  turbulent  patches  was  computed,  using  contours  of  the  finite¬ 
time  Lyapunov  exponent.  Edge  detection  techniques  of  the  map  yielded  the  propagation 
speeds  of  the  leading  and  trailing  edges,  which  agreed  with  the  literature.  In  addition, 
the  conical  transformation  was  verified  for  the  naturally-triggered  spot.  It  was  shown 
that  the  shape  of  the  various  patches  of  turbulence  can  collapse  into  a  self-similar  shape. 

Finally,  a  comparison  of  the  A2  vortex  identification  criterion  and  the  finite-time  Lya¬ 
punov  exponent  was  presented.  Agreement  in  the  prediction  of  both  techniques  was  only 
observed  near  the  edges  of  the  spots.  These  edge- vortices  play  a  significant  role  in  the 
global  dynamics  of  spot  spreading.  However,  the  internal  structure  of  the  spot,  namely 
the  hairpin  vortices,  were  only  captured  by  the  A2  criterion.  This  was  attributed  to  the 
relatively  short  time  scale  of  the  formation  and  breakdown  of  the  hairpin  vortices,  in 
comparison  to  the  backward  integration  time  used  in  the  computation  of  the  FTLE. 
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